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We investigate theoretically the optical properties of Nano- Wires (NWs) with cross sections having 
either discrete or cylindrical symmetry. The material forming the wire is birefringent, showing a 
different dielectric response in the plane and along the axis of the wire, which is typically the case for 
wires made of wurtzite materials, such as ZnO or GaN. We look for solutions of Maxwell's equations 
having the proper symmetry. The dispersions and the linewidths versus angle of incident light for 
the modes having high momentum in the cross-section plane, so called whispering gallery modes, are 
calculated. We put a special emphasis on the case of hexagonal cross sections. The energy positions 
of the modes for a set of azimuthal quantum numbers are shown. We demonstrate the dependence 
of the energy splitting between TE and TM modes versus birefringence. The polarization conversion 
from TE to TM with increase of the axial wave vector is discussed for both cylindrical and discrete 
symmetry. 

PACS numbers: 71.36.+c,71.35.Lk,03.75.Mn 



I. INTRODUCTION 

Nano- Wires are the objects with a cross-section dimen- 
sion a reduced to nanoscale values, several orders of mag- 
nitude smaller than their length L. The dynamics along 
the wire axis can be decoupled from the transverse one 
in case L » a and the system can be treated as quasi 
one-dimensional. Low dimensional systems in general re- 
veal novel phenomena and are excellent candidates for 
applications in new technologies. NWs could be used in 
communications^, quantum computation 2 -, or biological 
sensors 3 -. 

In the last years, significant progress has been achieved 
in the growth of semiconductor NWs. Many new inter- 
esting optical effects have been reported for such struc- 
tures, like ultraviolet lasing under optical pumping in 
ZnO NWs^ and recently, polaritonic effects. ID exciton- 
polaritons^ appear due to the strong coupling of exci- 
tons with photonic whispering gallery modes (WGMs)S. 
These modes owe their name to their similarity with 
acoustic resonances in real galleries. They propagate in 
the NW's cross-section and due to azimuthal momen- 
tum undergo multiple internal reflections. The number 
of these reflections can be very high, resulting in a high 
quality factor. 

The problem of resonant mode frequencies and of their 
life time in dielectric resonators has been studied for the 
cylindrical and hexagonal cross-sections both in meso- 
scopic (large ka)L and microscopic (small ka)^ regimes; 
where k is the wave vector of the incident light. NWs 
with cylindrical cross-section are very well theoretically 
studied in the isotropic case, being a textbook subject^. 
On the other hand, discrete cross-sections symmetries are 
much less studied whereas they are of particular inter- 
est because real structures often have polygonal cross- 
sections. Wurtzite or diamond crystals (like ZnO or 



GaN) generally form hexagonal NWs, but rectangular 
and triangular shaped NWs are also possible^. 

A remaining theoretical task in these discrete symme- 
try cross section systems is to calculate the dispersion 
of WGMs (dependence of their energy versus the axial 
wave vector), as well as their linewidth. The impact of 
the anisotropy of the dielectric response of wurtzite ma- 
terials on the polarization cigenstatcs of the WGMs has 
also not been addressed so far, to the best of our knowl- 
edge. On the other hand, it is well known that semi- 
conductors with wurtzite structure possess birefringent 
optical anisotropyii. NWs fabricated from such materi- 
als have directional dielectric response characterized by 
two refractive indexes - along the axis of anisotropy n z 
and perpendicular to it n r . Another source of birefrin- 
gence dn = n z — n r comes from axial variation of the NW 
radius due to inevitable growth imperfections. 

In this paper, we determine the optical eigenmodes 
of NWs of arbitrary cross-sections, considering both 
isotropic and anisotropic media. In the case of circu- 
lar NW cross section, the solutions are cylindrical har- 
monics of particular azimuthal quantum number m. For 
arbitrary discrete geometries, the appropriate linear com- 
bination of cylindrical harmonics is used to fulfill the 
boundary conditions for the tangential electric and mag- 
netic fields on the NW surface. The dispersions, line- 
widths, polarization, and field distribution of the modes 
versus the angle of incident light (longitudinal wave vec- 
tor) are deduced from this model. In case of birefringent 
media, the k z = WGMs are either TE or TM polarized. 
For other longitudinal wave vectors, strong mixing of po- 
larization occurs in both isotropic and anisotropic NWs 
leading to the formation of hybrid modes - EH and HE, 
depending on whether electric or magnetic field domi- 
nates in z-dircction. Nevertheless, in this paper we will 
use the notation TE for the former and TM for the later. 
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II. SYSTEM AND MODE SYMMETRIES 

We consider an infinitely long wire with its main axis 
along the z-direction and having either a circular or a 
polygonal cross section. The Figure [IJa) shows a finite 
length hexagonal cross section wire. The translational 
symmetry along the z-dircction imposes an additional 
factor in the wave functions having the phase form e tk * z . 
k z is a good quantum number which can take any real 
value for an infinitely long wire. 

The axial symmetry of regular polygons is the n-fold 
rotation C n , where n is the number of polygon sides. 
For a circular cross section n = oo, which corresponds 
to rotation by infinitcsimally small angles around the z- 
axis. The symmetry group Coo is continuous and, like in 
the case of translation, there is a good quantum number 
m G Z. Then, the tp- dependence of eigenmodes for each 
to is given by a single angular harmonic e mtp . 

For regular polygons, n is finite: triangle C3, square 
C4, pentagon C5, hexagon Cq etc. These systems re- 
main invariant for in-plane rotations by corresponding 
discrete angles 2irm/n. The existence of the mini- 
mal angle of rotation in these symmetries results in fi- 
nite number of irreducible representations of the group 
with to = 0, ±1, ±(n — 2)/2,n/2 (n even) or to = 
0, ±1, ±(n— 1) /2 (n odd)^£. The transformation prop- 
erties of higher to-s appear to be equivalent to those of 
the first Brillouin's zone in m-space. The eigenmodes be- 
havior along Lp coordinate is no more given by a single 
angular harmonic, but rather by their linear combination. 

Another important class of polygonal cross section 
symmetries are reflections - vertical mirror symmetries. 
The hexagon, for example, is unchanged under reflec- 
tions in the six vertical planes: 3 containing big hexagon 



diagonals and 3 connecting the centers of the opposite 
sides. These operations <7d,i and a v _i transform in-plane 
rotations by the angle tp,R(<p), into rotations by the op- 
posite angles -ip, R(—ip): <y v ^R{ l f>)<Jv,d = R(— y)- An 
important physical consequence is the equivalence of to 
and — to. The whole symmetry group describing the ax- 
ial symmetry of polygonal cross-section NWs containing 
both n-fold rotations and mirror reflections is denoted 
with C nv . 

The electric E = (E r ,E v ,E z ) and magnetic field 
H = (H r ,H V7 H z ) transform differently under re- 
flections. The electric field is a polar (ordinary) 
vector and magnetic field, being the cross-product 
H = (l//i)V x A, is a pseudovector. Reflec- 
tions change an eigenmode (E r , H r , E v , H v , E z , H z ) to 
(E r , —H r , —E^, H^, E z , —H z ) . It is interesting to see 
how the usual transverse electric (TE) and transverse 
magnetic (TM) eigenmodes for k z = are modified un- 
der reflection. TE modes transform from (E r , E v , H z ) 
to (Erj—Epj—Hz), i.e. their parity is opposite to that 
of the scalar function H z . TM modes transform from 
(H r , Hp, E z ) to (—H r ,H v ,E z ) after reflection, preserv- 
ing the parity of the scalar function E z . 

III. FORMALISM 

We start with the Maxwell's equations written in the 
frequency domain 



V x H(r,oj) = -i-D(r,u) 
c 



V x E(r,ui) = i-B(f,uj) 
c 



(1) 



(2) 
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FIG. 1: a) Hexagonal cylinder of finite height and the radius 
of its circum circle. The arrows show the plane of S and 
P polarized light. 6 is the angle of incidence b) Tangential 
component of the electric field on the surface of the hexagon 
at the point (r k ,4> k ) 



FIG. 2: Modes of a cylindrical cross-section nanowires versus 
angular momentum m (large symbols). The small symbols 
represent the modes of an hexagonal cross section nanowire 
which are coming form the superposition of the modes of the 
cylinder having an hexagonal symetry. 
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VD(r»=0 (3) 



VB(r» = (4) 

Taking the curl of the second equation and using 
B(r,uj) = fj,H(f,oj) from Eq.(p} follows the wave equa- 
tion 

uj 2 

V x V x E(r, lo) - -iruDff. uj) = 0. (5) 

& 

Displacement field is given by D(r,uj) = eE(f,iv), and 
we consider anisotropic permittivity e of the form e = 
diag(e r , s r ,s z ), where e r and e z are the permittivities in 
the NW cross section and along the z-axis respectively. 

The double curl of the electric field in the equation ([5]) 
can be rewritten as 

VxVx£ = V(VE) - AE. (6) 

In an isotropic medium this equation immediately 
gives VE = and Eq.© reduces to —AE. In an 
anisotropic medium, both terms of the equation © arc 
in general non-zero. However, one still can decompose 
the waves into pure TE ones (for which E z = and 
VE = kE = 0) and pure TM ones (for which H z = 
and VE = kE ^ 0). 

Let's now derive the equation for E z for the TM wave 
from equation (0. Using the equation ((3j> we calculate 
VE = dE z /dz(l — e z /e r ) which is non zero in case of 
TM waves in anisotropic materials. The z-component of 
equation J6]) reads 
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Inserting the equation (|7|) in the wave equation ([5]) 
allows obtaining the equation for the E z for the TM wave: 

The equation for the H z field for the TE wave has the 
form of the usual Hclmholtz's equation for an isotropic 
problem: 

AH Z + ^e r H z = 0. (9) 

So far we have considered infinite homogeneous me- 
dia. The solution for the wires can be obtained by 
matching the boundary conditions for the bulk waves 
on the wire surface. These boundary conditions mix 
bulk TE and TM waves, making decomposition into 
pure TE and TM modes only possible for k z = 0. Be- 
cause of the non-separability of TE and TM modes, we 
have to simultaneously solve equations © and ©. For 
cylindrical wires, we search for the solution in the form 
RE,H(r)e l ( mip+k * z ~ u ' t,> as it follows from the section on 
the mode symmetry. The radial parts Re,h(.t) satisfy 

!(r^) + (^-£)«„(r)-0, dO, 
where 

Qh = <te(£r/e z ) 1/2 - (e r c. 2 /c 2 - kl) 1 ' 2 (11) 

are in plane wave vectors, of magnetic and electric field 
within the wire. In the isotropic case qn = qE- Outside 
the NW, the same equation (|10p holds for the radial part 
with q out = (E ou t<jj 2 j<? — k 2 ) 1 / 2 ; e out is the dielectric 
constant of wire's environment that will be assumed to 
be isotropic. 

The equation (TIT))) is a Bessel's differential equation. 
The solutions inside the wire are the linear combina- 
tions of Bessel's functions of the first kind J m {x), whereas 
propagating solutions outside are the linear combinations 
of Hankel's functions H^x) of the first kind. Therefore 
one can look for the solution for E\ n and H™ fields in 
the following form: 

E\ n = ^A m J m (q E r)4> m -H^ 1 = ^ j B m J m {q H r)<t> m 

m m 

(12) 



FIG. 3: Cylindrical mode evolution versus the shape param- E° ut = V C m H^(q out r)(j) m ; H° ut = V D m H^{q out r) 
eter x which varies from for a cylinder to 1 for an hexagonal 



wire 



(13) 
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where <j) m = cxp(i(rmp + k z z — tot)). The transverse in- z-components (12,13) through the following equations: 
plane components of the fields can be deduced from the 

I 



q 2 E T = i-V T xH z + ik z \7 T E z , (14) 



q *H T = -i-V T xE z + ik z V T H z , (15) 
c 



where q = qn,E for internal fields and q = q ou t for external field. Inside the wire we directly find: 



E T = ~E ( A m 1 ^r 1 J m (qEr) +iB m — J m {q H r) ) <f> m , 
lH r c( 1h 



E™ = V ( iA m ^j' m (q E r) - B m -^J m (q H r) ) <f> m , (17) 

1h C lH r 



H?=J2( iA m ^j m (q E r) - B J^J m {q H r) ) m , (18) 



H l r n = V A m ^-J m {q E r) + iB m — J m {q H r) <j> m . (19) 
^ x cq H r c( 1h ' 



The outside fields read: 



E T = -Y / (c m ^H} n (q out r)+tD m —H^(q out r)) <b m , 

_ V liut r c 1out J 



H T = E (c m ^H^q out r) l D m ^Hl(q out r)) cj> m , 

^ \ c Qout lout 7 / 



(20) 



K Ut = E {c m —Bi (q out r) - il^-^-i&fo^r)') rn , (21) 

m \ lout c Qout r J 



(22) 



H° Ut = E ( C m ^p^J m {loutr) + iD m ^j m {q out r) ) cf> m . (23) 



I 

From these expressions, one can check that it is not of TE and TM modes is no more possible and the modes 

possible to have a pure TE mode, putting A m = (no are mixed. 

longitudinal electric field E = 0). For k z + 0, the mag- Tq calculate ^ mode digpersi one needs to im _ 

netic field is also present m the transverse plane and the , n , , ,., . c , , n , , . 

. ., .. , , _ , . . 1 pose the proper boundary conditions for both electric 

last two terms m the equations (HUH) are non-zero. Sim- and magnctic fidds rcsulting in the system of equa _ 

ilar argumentation holds for the existence of a pure TM ^ for ^ fidds amplitudes . Thc tangential ex _ 

mode. In the general case of k z + 0, a sharp separation ^ &nd electrk &nd nctic ficlds should 
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match on the wire boundary. If we denote F iniOUt = 
( E in,out^ E in,out^ j^.otrt jjpotrtyr boundary conditions 
are fulfilled if SF = Fi n — F out is zero on the surface 
of the wire. In the case of cylindrical wires of radius a, 
the fields for a single angular harmonics arc proper so- 
lutions, and the last statement is valid for each of them 



independently. Tangential components on the boundary 
are 93-field (equations [TS] and ITS)) and for a chosen to, 
on the boundary r = a we can write matrix equation 
F m (u>,k z ;a)X m = 0. The vector X m being the set of 
partial amplitudes X m = (A m , B mi C m , D m ) and matrix 
F m (ui, k z ; a) given by the following expression: 



F m {w,k z ;a) = 



Jm(qEa) 

H 



(<lEa) 





I 



cq h 



-H} n (q out a) 


-H}n(qouta) 





-H m {q uta) 
-H m (q out a) 
^- a Hl(q out a) 



cq 



(24) 



Dispersion relation oj m (k z ) is calculated from the con- 
dition that the determinant of a homogeneous system of 
equations must be zero in order to have a nontrivial solu- 
tion X rn 7^ 0. These dispersions can be found by efficient 
iterative numerical procedure^, based on the lineariza- 
tion of the matrix F m (u>, k z ; a) with respect to u>. 

The Figure [2] shows ( large simbols) the calculated 
mode energies at k z = for a cylindrical NW. For each 
angular momentum to, and for each polarization (E z , H z ) 
a series of modes characterized by their radial quantum 
number n r are showing up. For simplicity we show only 
the modes n r = and n r = 1 on the Figure [2] calculated 
for a radius a = 0.6/u.m and index of refraction n = 2.52. 

Before discussing the details of the calculation proce- 
dure for non-cylindrical NWs, it is helpful to consider 
qualitatively the expected modification of the mode spec- 
tra induced by the deviation of the NW shape from the 
cylindrical geometry. When NW cross section is just 
slightly deviated from a circle, the qualitative description 
of the modes can be obtained by perturbative considera- 
tions. The small deviation of the boundary cross-section 
results in the mixing between the cylindrical modes of 
different momenta. The physical reason of this mixing is 
the additional scattering of a cylindrical mode on a non- 
cylindrical boundary. The 'selection rules' of this scat- 
tering obviously depend on the symmetry of the NW. 
In case of an hexagonal perturbation e.g., that is invari- 
ant in respect to any rotation on multiples of 2n/6 , the 
changes of angular momentum are the multiples of 6. As 
a result, each cylindrical mode acquires a 'tail' of other 
harmonics whose amplitude depends on the energy mis- 
match between the energy of the main mode and that 
of the member of the 'tail'. The less is the energy mis- 
match, the stronger is the admixture of the other state 
to the initial cylindrical mode. The states with oppo- 
site to are degenerate in the absence of magnetic field 
and consequently can be strongly coupled by the surface 
perturbation, in case the angular momentum difference 
corresponds to that of boundary perturbation. 

To illustrate the mixing between the modes the Fig- 
ure [5] shows the structure of the cylindrical modes for 



small angular momenta with superposed harmonics with 
the periodicity of a hexagon. 

For the hexagon, a strong coupling of degenerate har- 
monics occurs for m— (—to) = 6/c, and results in a strong 
mixing of the modes with e.g. to = ±3 . As a re- 
sult, the to = ±3 doublet splits into a pair of singlet 
states. The dependence of this splitting on the hexago- 
nal perturbation is illustrated on the Figure [3] for a set 
of cross-sections smoothly varying from a cylinder (pa- 
rameters of [2]) to an hexagon with the same area of the 
cross-section. The radial position of the boundary for 
each angle was obtained by linear interpolation between 
the cylinder (a; = 0) and the hexagon (x = 1), where x 
is the shape parameter. 

Let us now discuss in more details the calculation pro- 
cedure for non-cylindrical NWs. We assume that the 
solutions inside and outside the wire can be written as 
linear combinations of the cylindrical harmonics. Bound- 
ary conditions in this case are more complicated. Single 
angular harmonics are no more proper solutions because 
they cannot match the tangential field boundary condi- 
tions on a boundary changing with the azimuth. The 
solution has to be written as an infinite combination of 
cylindrical harmonics as shown in the equations (| 16ti23|) . 
Moreover, because the normal to the boundary is not 
parallel to the radial vector, the tangential fields on the 
wire's surface contain both radial and azimuthal compo- 
nents ( Figure QJb)): 

Ei n '° ut (r k ,<p k ) = El n >° ut (r kl <p k )sm(Te r )e r + (25) 
E^ out {r k ,ip k )cos{fe v )e v . 

The matching of inside and outside fields should be 
realized at any angle < tp k <27r/n with corresponding 
radius r k = r(ip k ) along the wire boundary. The r and 
tp components in the above expression are given by the 
equations (|1TI21[) and f| 16I20[) . The magnetic field tan- 
gential to the NW surface is described in a similar way. 
In order to keep a convenient matrix description of the 
problem, we are going to consider boundary conditions 
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only on a finite number of points on the surface, and 
also to consider a finite number of cylindrical harmon- 
ics in the expressions of electric and magnetic fields. In 
this framework, boundary conditions can be expressed 
through the matrix equation F(oj,k z )X = 0, where the 
matrix F(u>,k z ) is given by: 



the angular momentum operator. Nevertheless, we will 
associate a number m corresponding to the angular mo- 
mentum of the principal harmonic to each mode, like in 
the case of a cylinder. The dispersion of these modes 
uj(k z ) can be found solving det(F) = or alternatively 
by a more efficient numerical procedu re 13 ' 14 . 



F(u,k z ) 



F u 



F kl 



F 



1 77 1 



F 



Fkr, 



XM 



FkM 



F 



F, 



Km 



... F 



KM 



(26) 



Here K is the number of points we take on the bound- 
ary and M is the number of harmonics we sum up. Ff~ m 
stands for an analog of F m (uj, k z ;r k ) - the matrices given 
by the expression (|24p with tangential (^-components re- 
placed by (|25|) . These sub- matrices describe contribu- 
tion of a single m-harmonic to the boundary condition 
at the surface point r k = r(tpk). The vector X is now a 
4M-dimcnsional colon (Xi, . . . , X m , . . . , Xm) t and X m 
is defined again like in the case of circular cross-section 
X m = (A m , B m ,C m , D m ). The meaning of the matrix 
([26]) is that, at any point r k on the NW surface, the 
same linear combination of cylindrical harmonics allows 
to verify the boundary conditions at these points. Taking 
the number of harmonics equal to the number of points 
on the boundary K = M allows to make the matrix (f2"rj)) 
square. In this case the eigenmodes ui(k z ) are found as 
the solutions of the system F(uj,k z )X = 0. Each eigen- 
mode of a n-polygonal system contains in addition to the 
principal harmonic m all other harmonics which add to 
to it by rule m + kn, where k is an non-zero integer. Such 
eigenmodes, resulting from summation over different har- 
monics, do not have a well defined angular momentum 
because of the fact that they are not the eigenstates of 



IV. RESULTS 

In this section we consider the important particular 
case of the hexagonal cross section and compare the 
eigenmodes with the ones obtained for a cylinder. The 
special case of the hexagonal cross section is of strong 
practical interest since it is realized experimentally by 
wires made of wurtzite semiconductors such as GaN and 
ZnO. The comparison with the cylindrical geometry is 
also particularly relevant since it is a much more sim- 
ple problem to solve. In practice, modeling of hexagonal 
NWs is often performed using a cylindrical description. 

We start by deducing empirically, how many cylindri- 
cal harmonics should taken into account in order to get 
a good precision on the energy of the eigenmodes. We 
consider a NW with a hexagonal cross section and a cir- 
cumcircle radius of 330 nm. For the FigurelH we consider 
an isotropic dielectric response with an optical index 2.5. 
We look for the eigen mode m = in the energy range 
3.05-3.25 cV versus the number of cylindrical harmon- 
ics M we take into account in the calculation. One can 
clearly see that the energy of the mode converges for 
large enough M. In all the following calculations, we 
will use M = 15. The energy of the modes is studied in 
the angular momentum space, the m-space (Figure E]). 
The modes of a hexagonal NW with a circumcircle ra- 
dius of 330 nm are analyzed for m = 1, ...12, and the 
energies m = 7, ...,12 arc shown with some small shift 
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FIG. 4: Convergence of the mode m = being in the energy 
range 3.05-3.25 eV versus the number of cyiindricai harmonics 
M used in the calculation 



FIG. 5: Energy position of hexagonal modes being in the 
chosen energy range versus their angular momentum. The 
circles show the modes belonging to the second "Brillouin 
zone" 
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with respect to the first six harmonics in order to ease 
the comparison. The numbers of modes and their en- 
ergies repeat themselves with a periodicity Am = 6 as 
expected from the discussion on mode symmetries for a 
hexagonal system. All relevant physical properties of a 
hexagonal NW can therefore be deduced looking into the 
first "Brillouin zone" placed between m = and m = 5. 
Next, we compare the energy, dispersion, and polariza- 
tion of hexagonal and circular NWs. Such comparison 
has been already performed in£, but only from the point 
of view of the energies of the eigenmodes at kz = 0. We 
are focused on the particular energies ( 3.1-3.7 eV, Fig- 
ure [6]) corresponding to the energy range in which semi- 
conductor's excitons couple to light, like for example in 
ZnO NWs£. In general, in order to directly compare the 
properties of the modes in circular and polygonal cross 
section NWs we have to consider structures having the 
same cross section area. If the radius of the cylinder is a, 
then the diagonal of corresponding n-sidc polygon can be 
found from the formula d = 2a(csc(2n/n)2w/n) 1 ^ 2 . The 
Figure E] shows the energy dispersions versus the angle of 
the incident light 9 = arcsin(fc z /fc) ( Figure (Ha)) of the 
modes of two NWs with circular and hexagonal cross- 
section respectively and having the same area. In both 
geometries the modes appear in polarization pairs TE 
and TM in a narrow energy range. The upper mode of 
the pair is characterized by an angular momentum m and 
is TM polarized, whereas the lower mode has an angular 
momentum m — 1 and is TE polarized. The modes of the 
hexagon appear at higher energies and have a slightly 
changed dispersion compared to the one of the cylin- 
der. By using slightly different index of refraction for 
the cylinder and the hexagon, it is possible to match the 
dispersion of one eigcnmode of the two different struc- 
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FIG. 6: The lower energy branch show the dispersion of TE 
and TM modes of circular cross-section wires: TE, m = 8 ( 
x symbols, blue) and TM, m = 9 ( + symbols, red) and the 
corresponding hexagonal modes: TE, m = 2 ( dashed blue 
line) and TM, m — 3 ( full red line). The same for the upper 
energy range replacing mbym + 1. 
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FIG. 7: Polarization degree of the lower energy modes of the 
Figure [5] hexagon modes m = 2 ( TE, red full line) and m — 3 
( TM, blue dashed line) and corresponding cylindrical modes. 
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FIG. 8: linewidths of the lower energy modes of the Figure[6] 
The index of refraction is n = 6.25 and the cross section area 
S = 0.09tt fim 2 



tures. This shows that in practice the eigenmodc of an 
hexagonal NW can be reasonably described by a simple 
model assuming a cylindrical geometry. 

Both in the cylinder and in the hexagonal wires, the 
modes are purely TE and TM only at k z =0. The Fig- 
ure [7] shows the dependence of the polarization degree 
P = {Ite — Itm)/{Ite + Itm) °f the modes versus the 
incidence angle. In the case of cylindrical wires p can be 
expressed directly through the coefficients of the exter- 
nal fields as Pm = (\D m \ 2 - \C m \ 2 )/(\C m \ 2 + \D m \ 2 ). The 
polarization degree of a TE mode decreases from 1 ap- 
proaching zero value at higher angles. The polarization 
of a TM mode changes very slowly, remaining close to — 1 
for all 8. The polarization mixing in hexagonal NWs has 
a bit different behavior in comparison with the circular 
geometry case. This difference is most obvious for TE 
modes. After being almost constant for a wide range of 
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FIG. 9: TE-TM splitting of hexagon modes at zero angle of 
incident light. The upper two modes of hexagonal NW are 
showmTE ( m = 3) and TM ( m = 4) 
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4 lower modes shown on the Figure [6j The linewidths of 
the hexagon modes arc of the order of 10 meV which 
is more than an order of magnitude larger than for the 
cylinder. It is even much larger at higher angles for 
which the linewidths of the modes of the hexagon re- 
main roughly constant whereas the ones of the cylinder 
drop by several orders of magnitude. This is due to the 
presence of corners in a hexagonal NW which arc respon- 
sible for higher losses resulting in larger linewidths than 
in a cylindrical wire. This was illustrated in reference^ by 
analysing the mode width dependence versus the round- 
ing of the corners of a hexagon. 

Another important feature of wurtzite materials is 
their optical birefringence. The optical index along the 
main c-axis, corresponding to the z-axis of the hexagonal 
NW differs from the one in the plane. The effect of the 
birefringence on the longitudinal-transverse splitting En 
(energy splitting of TE and TM polarizations) is shown 
on the Figure [9] for 6 = 0. Elt depends linearly on 5n 
and we show dependance on positive birefrigence like it 
is the case in ZnO wires^. Even a small birefringence 
leads to a significant splitting. It is therefore important 
to take the anisotropy into account in order to be able to 
reproduce realistic experimental situations. The radial 
dependence of the electromagnetic field density is shown 
on the Figure [TU] for the case of a hexagonal NW. 
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FIG. 10: Energy density radial dependence of the lower pair 
of hexagon modes:TE ( m = 3) (full red line) and TM ( m = 
4) (dashed blue line) 



angles, it starts to decrease significantly near 8 ps 30 and 
becomes even slightly TM polarized between 40 and 60 
degrees. On the other hand, the evolution of TM modes 
is similar to the case of a cylindrical structure. Such 
behavior results form the mixing of the hight-m wispcr- 
ing gallery harmonics with low-m ones with larger radial 
numbers. 

The main difference between the hexagon NW and the 
cylinder is the linewidth of the eigenmodes which is much 
larger for the hexagon. This is demonstrated on the fig- 
ure 7 which shows the mode linewidth versus angle of the 



V. CONCLUSIONS 

In conclusion, we have developed a method which al- 
lows to solve Maxwell's equations in NWs of discrete sym- 
metries, and even ones showing an anisotropic dielectric 
response. This method can be applied to any system 
having the cross section symmetry of regular polygons. 
It allows to find the eigenmodes of the structure (whisper- 
ing gallery modes) labeled by their angular, or pseudo- 
angular momentum in the case of non-cylindrical struc- 
tures. The dispersion (dependence of the energy on the 
wave vector along the wire axis) , polarization, linewidth, 
and radial densities of the modes are calculated for the 
cases of hexagonal and circular cross sections having the 
same area respectively. The modes in both cases appear 
to be quite similar, except from the point of view of the 
linewidth, which is much larger for hexagonal NWs. Wc 
have found some interesting polarization mixing effect 
with the transformation of TE modes close to k z = in 
modes mainly TM polarized. 

The authors would like to thank Maxime Richard, Le 
Si Dang, and Aurelien Trichet for fruitful discussions on 
this subject. 
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